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Abstract 

The three-loop contribution to the MS single-Higgs-doublet standard- 
model cross-section a(W^W£ — > ZlZl) at s = (5Mh) 2 is estimated 
via least-squares matching of the asymptotic Pade-approximant predic- 
tion of the next order term, a procedure that has been previously ap- 
plied to QCD corrections to correlation functions and decay amplitudes. 
In contrast to these prior applications, the expansion parameter for the 
W^W^ — > ZlZl process is the non-asymptotically-free quartic scalar- 
field coupling of the standard model, suggesting that the least-squares 
matching be performed over the "infrared" /i 2 < s region of the scale 
parameter. All three coefficients of logarithms within the three-loop term 
obtained by such matching are found to be within 6.6% relative error of 
their true values, as determined via renormalization-group methods. Sur- 
prisingly, almost identical results are obtained by performing the least 
squares matching over the /i 2 > s region. 

Within standard-model single-Higgs-doublet electroweak physics, the cross- 
section for the scattering of two longitudinal W's into two longitudinal Z's at 
very high energy takes the form 

87T 3 

<r[s, L{ji),g{ji)\ = ~g^H[s, L(», s(m)]> (1) 
where the scale-sensitive portion of the cross-section [1] , 
H[s,L{^),g^)] = 5 2 (/i){l+[-4i(/i)-10.0] g (//) 

+ [c 3 i 3 (fi) + c 2 L 2 ( f i)+c 1 L( f i)+ c ] g 3 (p) + ...}, (2) 
depends on the renormalization scale fi explicitly through the logarithm 

L(ji) = Intf/s) (3) 
and implicitly through the MS quartic scalar-field coupling 



12L 2 (» + 68.667LO) + ( 93.553 + \ln 



(4) 



The three-loop coefficients {co, ci, 02, C3} in (2) are presently unknown. The 
factor of Mh appearing explicitly in the two-loop term calculated in [1] is a 
scale-independent pole mass; the coefficients {co, Ci, C2, 03} can also exhibit de- 
pendence on this mass without acquiring additional /i-dependence. 

In this note, we utilise asymptotic Pade-approximant methods to predict 
these four coefficients. Such methods have already been applied to predicting 
next-order terms within the renormalization group functions of QCD [2,3,4], 
supersymmetric QCD [2,5], and massive scalar field theory [3,4,6], as well as 
next-order QCD corrections to scalar and vector current correlation functions 
[4,7] and various decay processes [8,9,10]. In all of these applications, the Pade- 
estimation procedures are tested (with surprising success) against either those 
higher-order terms already known from explicit calculation, such as renormal- 
ization group functions in scalar field theories [11,12], QCD [13], and supersym- 
metric QCD [14], or against those coefficients of logarithms [such as {01,02,03} 
in (2)] which can be extracted via renormalization group methods [7,8,9,10]. Be- 
low, we shall apply the latter testing procedure to estimates of the next-order 
terms {01,02,03} in the MS cross-section (1). 

Given a perturbative series of the form 1 + Rig + Rig 1 + Ry,g z + ... where 
i?3 is not known, as is the case in (2), asymptotic Pade-approximant methods 
can be employed to show that [4] 

Rf + R1R2 

This result is, of course, contingent upon the field-theoretical series exhibiting 
appropriate asymptotic behaviour. Its derivation (explicitly presented in [8]) 
follows from the 0{1/N) error anticipated from an \N\1\ Pade-approximant 
prediction of Rn+2 [15], a semi-empirical behaviour which is seen to characterise 
a number of field-theoretical applications even when N is small. 

For the case of the series (2), however, R\ is linear in L and R2 is quadratic 
in L. Consequently, the prediction (5) for R s corresponds to a rational func- 
tion of L incompatible with the degree-3 polynomial in L anticipated from (2). 
Clearly, a procedure is required by which predictions for the polynomial co- 
efficients {00,01,02,03} in (2) can be extracted from (5). In past applications 
where the same problem arises [7,8,9,10], one method employed is a least-squares 
matching of (5) to the form R3 = cq + ciL(/j,) + C2L 2 (p) + c^L 3 (fi) over the full 
perturbative domain of fi. For QCD calculations this domain is ultraviolet; e.g. 
in estimating three-loop QCD corrections to B — * X u l~T>i the matching is over 
the ultraviolet domain /i > m,b(p) [8]. For the expression (2), in which the 
perturbative expansion parameter is the non-asymptotically-free quartic scalar 
coupling \jjg(fj,), the appropriate domain for such a least-squares matching is 
infrared. 

Thus, to obtain predicted values for {co, 01,02,03} for a given choice of Mh, 
we choose a least-squares matching over the region fj, 2 < s, or alternatively 
< w < 1, where w = n 2 / s is the argument of the logarithm (3). From (2) and 
(5), this matching is achieved by minimizing the function 
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X (c ,ci,c 2 ,c 3 ) = 

j r 2i?3( w ) • 

dw -^3-7 — „ „ ; — ; - ; — r — Q) — c\ln{w) — ciln (w) — c 3 m (w) 
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(6) 



_Rf(w) + R 1 (w)R 2 {w) 
with respect to cq, c\, c%, C3, where R\{w) and are explicitly given in (2): 

fl 1 (io) = -4Zn(«;)-10.Q, (7) 

i? 2 H = 12Zn 2 (uj) + 68.667/n(w) + ^93.553 + j^ra h^)) ■ ( 8 ) 

The lower bound of integration w m i n in (6) would ordinarily be zero to encom- 
pass the full /1 2 < s range. However, we are compelled to consider a nonzero 
value of w m in in order to avoid any integrand poles, as discussed below. The 
expressions (1) and (2) are stated in ref. [1] to be accurate (within single-digit 
percent errors) only in the high-energy limit y/s > 5Mh- Although the projected 
linear and quadratic dependence of c\ and Co on ln{s/Mfj) could, in principle, 
be extracted via Pade methods ,[] the relatively small coefficient of this logarithm 
in the previous-order term (8) necessarily implies a similar insensitivity to this 
logarithm in Pade estimates of next-order terms. Consequently, we restrict our 
analysis here to the s = (5Mfj) 2 kinematic boundary of applicability for (1) 
and (2). With this choice, the integrand of (6) acquires singularities at 0.0552, 
0.0821, and 0.0896. Consequently, we choose w min — 0.09 to include virtually 
all of the integrable infrared region, and we find that 

X 2 (c ,ci,c 2 ,c 3 ) 

= 248033 + 813.779c - 332.772 Cl + 244.574c 2 - 242.664c 3 

+ 0.91c 2 - 1.38657c c! + 1.72946c c 2 

- 2.67527c c 3 + 0.864732c? - 2.67527cic 2 + 4.64965cic 3 

+ 2.32482e| - 8.67669c 2 c 3 + 8.48631^ (9) 

By then optimizing (9) with respect to cq, c±, c 2 , c 3 , we obtain the following Pade 
predictions for these coefficients: 

c Pade = _ 89 g 5 c Pade _ Q Pade = _ 2 g 8) c Pade _ _ 3Q 5 ^ 

As noted earlier, the true values of the coefficients c\, c 2 , c 3 can be extracted 
via the scale-invariance of the physical cross-section (1): 

O = ^^L{s,m,g(n)} = H + 0(g)^ (11) 



where [11] 



(3{g) = 2g 2 - ^g 3 + 27.803<? 4 + ... . (12) 



1 In ref. [10], the polynomial dependence of three- loop order terms in H — > gg on the 
logarithm of the pole-mass ratio Mjj/Mt is similarly extracted. 
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One can verify that the known terms in (2) satisfy the renormalization-group 
equation (11) to 0(g 3 ) and 0(g 4 ), as is evident from the series expansions 

OH 

— = -4 3 3 + (24L + 68.667).g 4 + (3c 3 L 2 + 2c 2 L + Cl )g 5 + ... , (13) 



/3( 5 )— = 4 ff 3 + (-24 J L- 68.667).g 4 + (96L 2 + 601.33L + 934.028 

+ 5.333;n(s/M|f)) 5 5 + ... . (14) 

We find upon incorporating 0(g 5 ) terms of (13) and (14) into the right-hand 
side of (11) that 



ci = -934.028- 5. min{s/M 2 H ) — ► -951.2, 

s=25M 2 H 

c 2 = -300.67, c 3 = -32. (15) 

The Pade predictions (10) for a, c 2 , c 3 are respectively seen to be within relative 
errors of 6.6%, 4.3%, and 4.7% of their true values (15). 

Curiously, the accuracy of these Pade results does not appear to be contin- 
gent upon the matching being performed over the "infrared" /1 2 < s = 25 Mfj 
range, as motivated by the non-asymptotically free character of the scalar field 
coupling g(/i). Indeed, a potential drawback of fitting over the fj, 2 < s (hence, 
w < 1) region, as in (6), is the negativity of ln(w) over the entire range of in- 
tegration. Since the Ci ultimately obtained in (10) are all same-sign (negative), 
cancellations necessarily occur between successive Ckln k (w) terms in the inte- 
grand of (6) in the best-fit region of Ck parameter-space. To address this issue, 
we have also performed a fit of the Pade-prediction (5) to the third-subleading 
order of (2) over the entire fi 2 > s (i.e. w > 1) region in which ln(w) is posi- 
tive. This entails integration of the integrand of (6) with appropriately modified 
bounds of integration to encompass the ultraviolet region: 



f di 



tr\ ..I 2 



dw[...} 2 . (16) 
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We then find that 



X 2 (c ,ci,c 2 ,c 3 ) = 1.45665 • 10 7 + 5095.28c + 10294.6ci 



+ 35519.6c 2 + 167166c 3 + c\ + 2c oCl 



+ 4c c 2 + 12c c 3 + 2c 2 + 12cic 2 + 48cic 3 
+ 24c 2 , + 240c 2 c 3 + 7204 (17) 

which upon optimization, yields values for Cfe[c = —896, c\ = — 889, c 2 = 
—289 c 3 = —30.9] that are virtually the same as those listed in (10). The 
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small relative errors characterising our Pade estimates of Ci, c 2 , C3 are compara- 
ble to those characterising Pade estimates of renormalization-group accessible 
coefficients within next-order QCD corrections to other processes [7,8,9,10], and 
suggest similar accuracy in the estimated value of the renormalization-group 
inaccessible three-loop coefficient c in (10). 

We therefore conclude that Pade-approximant predictions of the next order 
contribution to WW — > ZZ at very high energies appear to be consistent and 
reliable. It should also be noted that higher order /3-function terms associated 
with the evolution of the quartic scalar-field coupling constant (4) are them- 
selves accurately predicted by the same asymptotic Pade-approximant methods 
employed above for <r(WW — > ZZ). If we express the /3-function (12) in the 
form 

13(g) = 2g 2 (l + R l9 + R 2 g 2 + R 3 g 3 + ...); Rt = -13/6, R 2 = 13.915, (18) 

we predict via (5) that R3 = —133.6, or alternatively, that the predicted next 
term in the series (12) is — 267. 2g 5 . This is quite close to — 266.495g 5 , the 
true calculated value [11] of the next-order /3-function contribution. Similarly 
close agreement between the somewhat more complicated asymptotic Pade- 
approximant prediction and the explicit calculation of the 0(g 6 ) contribution 
to this /3-function is demonstrated in ref. [6]. 

Thus, the results presented above are an example of how Pade estimation 
procedures can anticipate next-order contributions whose exact values are ob- 
tainable only by lengthy calculation. For the particular process in question, the 
distinction between two- and three-loop order results is seen to be unimportant 
unless the mass of the (Salam- Weinberg) Higgs field mediating the scattering 
process is very large. This is illustrated in Figures 1-3, which compare two loop 
and three loop expressions for the scale sensitive portion (2) of the WW — > ZZ 
cross section (1) for Higgs-field masses of 200, 400, and 600 GeV, respectively. 
Only for the largest of these three choices is an appreciable difference anticipated 
between two- and three-loop order predictions for the cross section. 
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Figure 1: The two-loop (bottom curve) and predicted three-loop (top curve) 
expressions for the scale-sensitive portion H[s — (5Mff) 2 ,L(/i),g(/j)] of the cross- 
section (1) are plotted for Higgs mass Mh = 200 GeV. 
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